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Abstract 

Background: Arsenic (As) is a toxic metalloid found ubiquitously in the environment and widely considered an 
acute poison and carcinogen. However, the molecular mechanisms of the plant response to As and ensuing 
tolerance have not been extensively characterized. Here, we report on transcriptional changes with As treatment in 
two Arabidopsis accessions, Col-0 and Ws-2. 

Results: The root elongation rate was greater for Col-0 than Ws-2 with As exposure. Accumulation of As was lower 
in the more tolerant accession Col-0 than in Ws-2. We compared the effect of As exposure on genome-wide gene 
expression in the two accessions by comparative microarray assay. The genes related to heat response and 
oxidative stresses were common to both accessions, which indicates conserved As stress-associated responses for 
the two accessions. Most of the specific response genes encoded heat shock proteins, heat shock factors, ubiquitin 
and aquaporin transporters. Genes coding for ethylene-signalling components were enriched in As-tolerant Col-0 
with As exposure. A tolerance-associated gene candidate encoding Leucine-Rich Repeat receptor-like kinase VIII 
(LRR-RLK VIII) was selected for functional characterization. Genetic loss-of-function analysis of the LRR-RLK Vlii gene 
revealed altered As sensitivity and the metal accumulation in roots. 

Conclusions: Thus, ethylene-related pathways, maintenance of protein structure and LRR-RLK Vlll-mediated 
signalling may be important mechanisms for toxicity and tolerance to As in the species. Here, we provide a 
comprehensive survey of global transcriptional regulation for As and identify stress- and tolerance-associated genes 
responding to As. 
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Background 

Arsenic (As) is a ubiquitously present non-essential metal- 
loid of serious environmental concern because of ever- 
increasing contamination [1,2]. Naturally high levels of As 
in drinking water have caused major human health prob- 
lems in the United States, Argentina, Taiwan, and most 
notably Bangladesh and India, where tens of mUlions of 
people have been affected [3,4]. As finds its way into the 
food chain through irrigation with contaminated ground- 
water [5]. Uptalce of As in plant tissues adversely affects 
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the plant metabolism and leads to a significant reduction 
in crop yield [6-8]. Understanding of As-induced stress 
and ensuing tolerance would be beneficial for the develop- 
ment of As-resistance crops and other economically im- 
portant plants [8]. 

Roots are involved in mineral acquisition by plants, and 
function at the interface with the rhizosphere. Alterations 
of root architecture and inhibition of root elongation 
are considered primary symptoms of As-toxicity [9,10]. In 
many circumstances, it is the As -sensitivity of the root 
that limits the productivity of the entire plant [11]. Hence, 
plants exposed to As show inhibited root growth and re- 
duced photosynthetic rate [8,12]. While plant roots are 
the first organs in contact with As, assaying the processes 
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occurring in the roots could provide potential strategies to 
determine how plants respond and adapt to the heavy 
metal stress. 

As is present both as As (III) and As (V) in the environ- 
ment, with As (V) being more prevalent than As (III) in 
soils [8,13]. The mechanism by which As is taken up by 
plants differs for As(III) and As(V) [14]. As (III) uptal<e is 
thought to occur through the aquaporins of roots [14,15]. 
Higher plants take up As (V) as the dominant form of 
phytoavaOable As in aerobic soils. As (V) is a phosphate 
analogue and easily incorporated into plant cells through 
the high-affinity phosphate transport system [8,14]. As (V) 
replaces phosphate in ATP to form unstable "ADP-As" 
complexes and disrupts energy flow in cells [8]. As (III) re- 
acts with the sulfhydryl groups of enzymes and proteins, 
thereby inhibiting cellular function and resulting in cell 
death [8,16]. In addition. As stimulates the formation of 
free radicals and reactive oxygen species, thus resulting in 
oxidative stress [9,10,17]. 

Detoxification and tolerance mechanisms often involve 
extrusion of the toxic ions from cells, sequestration within 
internal organelles, reduced toxin uptake, and chelation by 
metal-binding proteins such as phytochelatins (PCs) [18]. 
These mechanisms reduce the cellular content of the toxic 
agent, although the molecular basis may differ among 
metals and organisms [19-21]. Differences in As sensitiv- 
ities exist among plant accessions or varieties [11,22,23]. 
Understanding the mechanisms underlying reduced As 
sensitivity and genes responsible is important in directing 
future breeding to counter As stress. The potential of As 
tolerance, based on such a concerted response of the vari- 
ous pathways, would also depend on an early perception 
of As-induced stress [10]. The precise mechanisms of 
perception of As-induced stress in plants remain to be 
elucidated. 

Reducing As intake requires identifying the mecha- 
nisms implicated in As toxicity and tolerance in plants. 
A substantial number of genes are differentially regu- 
lated by As stress in various plant species [9,10,24]. 
However, genes identified as As responsive are in part 
related to a general stress response resulting from the 
toxic effects of As and are unlikely to play a significant 
role in As tolerance. High-throughput gene expression 
profiles with microarray technology and their application 
in comparative studies help to reveal the roles of differ- 
ential gene regulation in As toxicity and tolerance. For 
more insight into the molecular basis of As toxicity and 
tolerance responses in plants, we performed transcrip- 
tional profiling of 2 contrasting accessions of Arabidop- 
sis thaliana, Col-0 and Ws-2. We focused on genes that 
are commonly and specifically regulated by the two ac- 
cessions and discuss the putative functions of identified 
genes in several biochemical pathways in terms of As 
toxicity and tolerance. New genes identified may provide 



more clues to understanding the molecular mechanism 
of response to As-induced stress in plants. 

Results 

Effects of As stress on root elongation in Arabidopsis 
accessions Col-0 and Ws-2 

Root growth inhibition is the primary response of the plant 
exposed to heavy metals. We analyzed the effect of As (V) 
exposure on primary root elongation to evaluate As tox- 
icity and tolerance in A. thaliana accessions, including 
Col-0, Ws-2 and Ler-0. As a result, we examined acces- 
sions for As tolerance and exposed 1 tolerant (Col-0) and 
1 sensitive (Ws-2) accession to As for 2 d for measuring 
root elongation (Figure la). At 100 |iM, As significantly re- 
duced root elongation in both Arabidopsis accessions as 
compared with normal conditions (Figure la). Of note, 
Col-0 showed less reduction in root elongation (15%) than 
Ws-2 (50%). The root elongation decreased with increas- 
ing As concentration (Figure lb). At 200 |iM As, the re- 
duction in root elongation in Col-0 was 40% as compared 
with 60% in Ws-2. At 300 ^iM As, root growth was min- 
imal in Col-0 and completely inhibited in Ws-2. As at 200 
and 100 \iM. inhibited root elongation by approximately 
one-half in Col-0 and Ws-2, respectively. Col-0 was more 
tolerant to As than Ws-2. Therefore, we hereafter refer to 
Col-0 as As tolerant and Ws-2 as As sensitive. 

Accumulation of As concentration in Arabidopsis roots 
with As exposure 

To investigate whether the differences in As tolerance in 
the two accessions was associated with As concentration 
in roots, we exposed root tips to As concentrations for 3 h 
and 48 h (Figure Ic). When exposed to As for 3 h. As was 
increased in roots of both accessions with increasing As 
concentration. The concentration of As was lower in the 
As-tolerant than As-sensitive accession (Figure Ic). Note 
that at 48 h, concentration of As was accumulated to a 
similar level in As-sensitive and As-tolerant accession, 
which possibly resulted from a more severely inhibited 
transpiration rate. Taken together, As accumulation in the 
roots of both accessions was associated with As concen- 
tration, with greater accumulation in Ws-2 than Col-0. 
The data may explain the differences in As toxicity and 
tolerance between the two accessions. 

Global expression profiles of Col-0 and Ws-2 in response 
to As 

To identify genes associated with As toxicity and tolerance 
in Arabidopsis, we used large-scale expression profiling. 
The As exposure data (Figure 1) helped determine a suit- 
able As concentration for microarray analysis. Treatment 
of Col-0 and Ws-2 with As at 200 and 100 ^M, respect- 
ively, had similar effects on root growth inhibition, indi- 
cating equal toxic strength. RNA samples were collected 
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Figure 1 Effect of As (V) treatment on the root elongation and metal accumulation of Arabidopsis thaliana accessions Col-0 and Ws-2. 

(a) Seedlings with different accession bad<grounds (Col-0 and Ws-2) were grown on quarter-strengtin MS medium for 4 d and tlien transferred to 
medium with 100, 200 and 300 (iM As (V) and grown for an additional 2 d (Scale bar, 2 cm ). As tolerance was determined by relative root growth 
after treatment, (b) Root length of plants was measured after treatment with As. Root samples were collected from 3 independent experiments 
(each from a pool of 7 root samples). Data are mean±SD. *P < 0.05 compared with Col-0 in each concentration of As. The difference in root 
elongation is significant according to Student's t test, (c) Accumulation of As (V) by Arabidopsis accessions Col-0 and Ws-2 was analyzed with 
ICP-AES. The 4-d-old seedlings with different accession backgrounds (Col-0 and Ws-2) were transferred to medium containing 100 and 200 pM As 
(V) for 3h and 48 h, then root tips were collected and measured. Data are mean±SD calculated from 3 biological replicates per treatment. 
*P < 0.05 compared with Col-0. 



from root tips early (1.5 to 3 h) after As treatment to 
examine rapid changes in global patterns of gene expres- 
sion. We pooled RNA isolated from with 1.5- and 3-h As 
treatment to maximize gene discovery. 

Fold-change values were compared with a control sam- 
ple without As treatment. Differentially expressed genes 
were defined as those with at least 2-fold change in tran- 
script abundance and with an adjusted P value less than 



0.05 [25]. In general, the As-tolerant and As-sensitive ac- 
cessions showed similar expression patterns in the control 
conditions. Therefore, constitutive gene expression was 
similar between the two accessions before As treatment. 
In total, 620 probe sets were specifically upregulated in 
As-tolerant Col-0 plants exposed to 200 As (Col-0 
200 i^M As versus Col-0 control) and 756 were downregu- 
lated (Figure 2a). These 1376 As-regulated genes were 
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Figure 2 (See legend on next page.) 
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(See figure on previous page.) 

Figure 2 Microarray expression of genes in Arabidopsis Col-0 and Ws-2 plants exposed to As stress, (a) Venn diagram of regulated 
As-responsive genes extracted by comparing microarray probe sets of the 2 Arabidopsis accessions. The number of overlapping and non-overlapping 
genes early (1.5 to 3 h) after treatment with 100 or 200 |jM As is shown. The probe sets were seleaed on the basis of an adjusted P-value of <0.05 and 
a >2-fold change in gene expression. The area of the diagram is proportional to the number of genes that are up- or downregulated in response to As 
stress. To clearly differentiate As-regulated genes in Col-0 from that in Ws-2, the Col-O-specific set contains only those genes with both >2-fold change 
in abundance (compared with control treatment) in Col-0 and <1 .4-fold change in Ws-2. Likewise, the Ws-2-specific set contains genes with both > 2- 
fold change in abundance in Ws-2 and < 1 .4-fold change in Col-0. The general As-reguiated gene set contains genes displaying > 2-fold change in 
abundance in Col-0 or Ws-2 and >1 .4-fold change in another accession, (b) Displayed are genes associated with ubiquitin pathways and abiotic stress 
responses using MapMan software. Both sets of material were harvested from roots tissues treated with or without As stress (100 and 200 pM for Ws-2 
and Col-0, respectively). Red and blue signals represent a decrease and increase in transcript abundance, respectively, relative to water-treated control 
samples. The scale used for coloration of signals (log2 ratios) is presented, (c) Validation of representative As-induced genes by semi-quantitative 
RT-PCR analysis. Total RNA was extracted from root tissues of Arabidopsis plants with different accessions with As treatment. The As-treated roots were 
harvested at 1 and 3 h. The samples were pooled together. The transcript level of actin served as an equal loading control. 



unique to Col-0. In Ws-2, 59 probes were upregulated 
in response to 100 |iM As (Ws-2 100 ^iM As versus 
Ws-2 control) and 15 were downregulated. These 74 As- 
regulated genes were unique to Ws-2. A total of 558 probes 
were significantly upregulated in both Col-0 and Ws-2 
and 517 were downregulated in both (Additional file 1: 
Table SI for full list). Thus, a relatively larger number of 
genes showed changed expression with As in Col-0, 
which suggests differential cellular response with early 
As exposure. 



GO analysis of As-responsive genes 

We subjected genes up- and downregulated with As stress 
to GO analysis using the EasyGO program to understand 
the metabolic and regulatory differences between the two 
accessions on As exposure. The major GO terms for bio- 
logical processes for As-responsive genes are summarized 
(Table 1). Upregulated genes for both accessions were in 
functional categories corresponding to responses to heat, 
oxidative stress, and metal ion and carbohydrate stimulus. 
Col-0-specific genes were involved in responses to ethylene 



Table 1 Gene ontology categories corresponding to genes differentially regulated by As in Arabidopsis accessions 
exposed to As 



Biological process 



Accession 



No. genes^ 



Background 



False discovery rates 



Upregulated generally by As (V) 

Heat response 

Oxidative stress 

Chitin response 

Metal ion response 

Carbohydrate stimulus 

Upregulated specifically by As (V) 

Ethylene response 

ABA response 

Heat acclimation 

Osmotic stress 

Toxin response 

Downregulated generally by As (V) 

Cell wall organization 
Secondary metabolic process 
Cytokinin response 
Glycoside biosynthesis 
Transport 

Downregulated specifically by As (V) 

Receptor protein signaling 



Col-0, Ws-2 
Col-0, Ws-2 
Col-0, Ws-2 
Col-0, Ws-2 
Col-0, Ws-2 

Col-0 (Ws-2) 
Col-0 (Ws-2) 
Col-0 (Ws-2) 
Ws-2 (Col-0) 
Ws-2 (Col-0) 

Col-0, Ws-2 
Col-0, Ws-2 
Col-0, Ws-2 
Col-0, Ws-2 
Col-0, Ws-2 

Col-0 



49,32 
42,29 
41, 11 
44, 17 
50, 14 

16,(5) 
14,(6) 
5(2) 
16(7) 
16 ( 10 ) 

49, 12 
47, 18 
8,6 
12,11 
73,33 

30 



114 

258 
127 
367 
203 

143 

262 

15 

391 

66 

263 

359 

75 

32 

1732 

130 



1.70E-27, 1.95E-25 
1.06E-9, 2.87E-13 
7.05E-19, l,39E-5 
2.58E-06, 0.00799 
5.80E-21, 0.000626 

5.01 E-05 

0.00773 

0.00154 

0.0346 

5.54E-12 

3.01 E-1 2, 3.81 E-05 
2.50E-07, 9.42E-06 
0.0196, 0.0066 
8.31 E-06, 4.30E-11 
0.0499, 0.0289 

l,96E-09 



"Accession specificity were identified according to false discovery rates (FDR < 0.05) of gene ontology analysis. Accession and no. of genes in parentheses refers to 
FDR > 0.05 in each gene ontology term. 
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Stimulus and abscisic acid (ABA) and heat acclimation. 
Ws-2-specific genes were involved in osmotic stress and 
toxin response. Downregulated genes in both accessions 
were predominately involved in cell wall organization, re- 
sponse to cytokinin, glycoside biosynthesis and transport 
(Table 1). Col-O-specific downregulated genes were related 
to receptor protein signaling pathway with As stress. 

To expand the functional significance of the As- 
responsive genes in Col-0 and Ws-2, we used MapMan 
representations to highlight microarray data and bio- 
chemical pathways. Col-O-specific genes with changed 
expression were involved in the proteasome and heat 
stress (Figure 2b). The sulfur assimilation pathway, which 
leads to phytochelatin biosynthesis, is affected by As 
stress in both accessions (Additional file 2: Figure SI). 
The gene encoding phytochelatin synthase was upregu- 
lated only in Col-0 plants after exposure to As stress. 

Oxidative stress-related genes regulated by As exposure 

The As-responsive genes can be grouped into various 
biological processes such as oxidative stress, transporter, 
hormone homeostasis and signal transduction. First, the 
expression levels of the genes dealing with the responses 
to oxidative stress were compared between the two ac- 
cessions (Additional file 1: Tables S2 and S3). Expression 
of genes coding for alternative oxidase, dehydroascor- 
bate reductase, glutaredoxin, peroxidase were regulated 
in both accessions on exposure to As stress (Additional 
file 1: Table S2). In addition, a number of genes encod- 
ing thioredoxin were mostly induced in tolerant Col-0. 
Genes encoding glutathione S -transferases (GSTs) were 
predominately upregulated by As in both accessions. A 
total of 14 genes encoding for GST were upregulated in 
As-sensitive Ws-2, and 9 GST genes were induced in 
As-tolerant Col-0. Most of the As-responsive GST genes 
belong to the 40-member Tau class GSTs [26]. Fewer 
GST-related genes were upregulated in As-tolerant Col-0 
than in As-sensitive Ws-2. 

Transporter genes regulated by As exposure 

A significant number of genes encoding transporters were 
differentially expressed in the As-sensitive and As-tolerant 
accessions with As exposure (Additional file 1: Table S4). 
ATP-binding cassette (ABC) transporters comprised the 
largest group of transporter-related genes (Additional file 1: 
Table S2). A total of 9 and 14 ABC transporter genes were 
differentially regulated in the Ws-2 and Col-0 accessions, 
respectively. Transcripts for genes from the multidrug and 
toxic compound extrusion (MATE) transporters and anti- 
porters were also regulated by As treatment. Genes an- 
notated as aquaporin, LeOPTl oligopeptide transporters 
(OPT) and sugar transporter were predominately down- 
regulated in Col-0 with As exposure, which suggests the 



differential regulation of transporter-related genes between 
the two accessions in response to As stress. 

Hormonal genes regulated by As exposure 

To understand the expression pattern of genes involved 
in hormone pathways between the As-sensitive and 
As-tolerant accessions, we analyzed transcripts related to 
hormone metabolism and found a potential role of ABA, 
brassinosteroid, cytokinin and ethylene in an early re- 
sponse to As stress (Additional file 1: Table S5). Genes en- 
coding proteins involved in cytokinin homeostasis were 
mostly downregulated by As stress in both accessions 
(Additional file 1: Table S2). ABA-related genes were regu- 
lated (expression level 2 ~ 5 fold change) predominately in 
As-tolerant Col-0. In total, 14 and 5 ethylene-related genes 
(expression level 2 ~ 26 fold change) were regulated in 
the Col-0 and Ws-2 accessions, respectively. Ethylene 
biosynthesis- and signaling-related genes were regulated 
by As stress, especially in As-tolerant Col-0. 

Transcription factors and protein kinases regulated by As 
exposure 

In total, 200 and 69 genes encoding transcription fac- 
tors were differentially expressed with As exposure in 
As-tolerant Col-0 and As-sensitive Ws-2, respectively 
(Additional file 1: Table S6). Col-0 showed more genes 
induced: APETALA2I ethylene-responsive-element-binding 
protein {AP2IEREBP), AuxlIAA, heat shock transcription 
factor {HSF), and WRKY (Additional file 1: Table S2). 
Genes encoding bHLH, C2H2, GARP-G2-like, and MYB 
were predominately downregulated in Col-0. Further, 
members of the AP2/EREBR HSF and WRKY families 
were upregulated specifically in Col-0 with As stress. 
Kinases may act as signals on the transcription factors, 
leading to the production of stress proteins and secondary 
metabolites that can act as either damage-causing or 
stress-countering agents. Genes involved in mitogen- 
activated protein kinase (MARK) and SNFl-related 
kinases (SnRKs) and leucine-rich repeat receptor-like 
kinase VIII (LRR-RLK VIII) pathways were more regu- 
lated in Col-0 than in Ws-2 by As stress (Additional file 
1: Tables S2 and S7). 

Identification of putative candidate genes for As 
tolerance in Arabidopsis 

We aimed to isolate the Arabidopsis genes responsible for 
As tolerance. The microarray data corresponding to the 2 
contrasting accessions and the dose-response effect were 
integrated to identify tolerance-associated genes. Accord- 
ingly, gene-fQtering criteria were based on 2-fold change in 
gene expression in As-treated Col-0 versus Ws-2 (Col-0 
200 [iM As/Col-0 Control versus Ws-2 100 ^iM As/Ws-2 
Control). Besides, a minimum 2-fold change in the dose- 
response effect in Col-0 (200 |iM As/Control versus 
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100 |iM As/Control) was included in the filtering cri- 
teria. We found 63 putative genes grouped according to 
biological processes by GO analysis (Table 2). These 
genes were further classified into regulatory genes, such 
as those encoding proteins responsible for signal trans- 
duction, transcriptional regulation, GTP binding and 
the proteasome-related pathway, and functional genes 
contributing to responses to heat, mitochondria elec- 
tron transport and responses to oxidative, salt and other 
stresses (Table 2). We found 12 putative regulatory genes 
encoding the proteasome-related pathway associated with 
As tolerance (Table 2). These genes belong to various 
types of proteasome-related components such as RING/ 
U-box superfamily protein (AT1G14200), RING-H2 
finger protein RHAla (AT4G11370), ubiquitin fusion deg- 
radation UFDl family protein (AT4G1542) and mem- 
brane-anchored ubiquitin-fold protein 4 (AT3G26980). In 
addition, tolerance to As was associated with the increased 
expression of genes related to heat stress, such as the 
heat-shock-protein 20-like chaperone superfamily protein 
(AT1G52560) and //S'f£2A (AT5G62020). 

A set of the tolerance-associated genes such as HSF 2A 
(AT2G26150), HSF B2A (AT5G62020), ethylene response 
factor (AT2G33710), GT-like trihelix DNA-binding protein 
(AT1G76880), ATSDG37 (AT2G17900) and LRR-RLK 
VIII (AT3G09010) were selected for validation by semi- 
quantitative reverse-transcription polymerase chain reac- 
tion (RT-PCR). The expression level of these genes was 
increased with As concentration in Col-0 (Figure 2c). In 
addition, the induction of the genes by As treatment was 
higher in Col-0 than in Ws-2. The regulation of these 
genes in response to As stress demonstrated the does- 
response and accession-specific effects. These results were 
consistent with the microarray data. 

Transcriptional characteristics of putative As-tolerance 
associated genes encoding LRR-RLK VIII 

Most of the As-tolerance associated genes were involved 
in signal transduction and regulatory mechanisms (47 
out of 63, Table 2). Their up-regulation at the early stage 
endorsed the trigger of downstream components to cope 
with the stressful condition. These regulatory genes may 
act as the fate dominators of As tolerance. Very little is 
known about the role(s) of LRR-RLK VIII gene in regula- 
tion of plants responding to heavy metal stresses. To 
characterize an As tolerance-associated gene candidate 
encoding LRR-RLK VIII, we analyzed transcriptional 
regulation for Arabidopsis LRR-RLK VIII genes regarding 
its response to As stress. According to the microarray 
data, four genes (AT1G53430, AT1G53440, AT3G09010 
and AT5G01950) belonging to the LRR-RLK VIII family 
were significantly upregulated by As stress (Table 2 and 
Additional file 1: Table S7). The data was further vali- 
dated by semi-quantitative (Figure 3a). The expression 



level of the LRR-RLK VIII genes (AT1G53430, AT1G53440 
and AT3G09010) was induced more significantly in As- 
treated Col-0 plants as compared to Ws-2. Next, we 
examined the expression of the LRR-RLK VIII genes 
(AT1G53430, AT1G53440 and AT3G09010) in response 
to various stresses (Figure 3b). The gene expression was 
strongly induced in root tissues with treatment of 
200 [iM As in comparison with H2O2 (100 ^iM) and Cu 
(25 |aM) stresses (Figure 3b). Therefore, the increase in 
LRR-RLK VIII gene expression was regulated in root tis- 
sues by As. The data demonstrated the specificity of the 
LRR-RLK VIII gene expression in response to As stress. 

Functional analysis of the LRR-RLK VIII gene in response 
to As stress using Arabidopsis T-DNA mutants 

To examine the function of an As tolerance-associated 
gene candidate encoding LRR-RLK VIII, we analyzed T- 
DNA mutant lines for Arabidopsis LRR-RLK VIII genes re- 
garding its response to As stress. Arabidopsis mutants 
with a T-DNA insertion in the locus coding for LRR-RLK 
VIII was characterized and subjected to As treatment 
(Figure 4a). The insertion of the mutant line (AT1G53440: 
SALK_057812) located 403 bp upstream of the start 
codon. The insertion abolished As-induced LRR-RLK VIII 
gene expression (Figure 4a). Two of the mutant lines for 
AT1G53440 (SALK_057812 and SALK_148231) showed 
alternation in As sensitivity (Figure 4b). At 200 (iM, As 
significantly inhibited root elongation in both the wild- 
type plants and the LRR-RLK VIII mutants when com- 
pared to untreated controls (Figure 4b). However, the root 
growth inhibition was less in the LRR-RLK VIII mutant 
lines as compared to the wild-type plants. Thus, knock- 
out of LRR-RLK VIII gene caused decrease in As sensitiv- 
ity. The level of As concentration in roots was less in the 
LRR-RLK VIII mutant line (AT1G53440; SALK_057812) 
than wild- type plants (Figure 4c). In addition, the func- 
tional specificity of the LRR-RLK VIII gene was assessed 
by treatment of the Arabidopsis mutant lines with various 
stresses such as As, Cu, Cd, Zn, H2O2 and salt stress 
(Figure 4d and Additional file 3: Figure S2). Cu, Cd, Zn, 
and salt stresses inhibited root growth to a similar level 
between wild-type plants and the LRR-RLK VIII mutant 
line (AT1G53440: SALK_057812). Alleviation of As- and 
H2O2- induced root growth inhibition was observed in 
the LRR-RLK VIII mutant line as compared to wild-type 
plants (Figure 4d). Thus, the loss-of-function mutation of 
LRR-RLK VIII caused As hyposensitivity, indicating its 
functional significance and specificity in response to As 
stress. 

Discussion 

Transcriptome profiling of As toxicity and tolerance 

We aimed to compare changes at the physiological and 
global gene expression levels of Arabidopsis accessions 
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Table 2 Selection of putative As tolerance-associated genes by comparing the expression ratio between the two 



accessions 

AGI Product encoded Tolerance ratio^ Dosage ratio'' 

Col-0 200/Ws-IOO Col-0 200/100 ((JIVI) 

Transcription factor 

AT3G17611 RHOMBOID-like protein 14 (RBL14) 2.84 2.12 

AT1G76880 GT-lil<e trihelix DNA-binding protein 6.41 2.11 

AT2G17900 ATSDG37 3.39 2.50 

AT5G62020 Heat 5hocl< transcription factor B2A 4.49 341 

AT2G33710 ERF (etliyiene response factor) transcription factor family 3.98 2.15 
Responses to heat 

AT5G59720 ATHSP18.2 30.66 2.18 

AT1G52560 HSP20-like chaperones superfamily protein 18.23 3.16 

AT5G37670 HSP20-like chaperones superfamily protein 12.70 3.06 

AT4G10250 ATHSP22 8.58 3.84 

AT5G15450 HSPlOO/ClpB 8.12 2.42 

AT2G32120 ATHSP70-2 4.51 2.89 

AT4G25200 AtHSP23.6 4.33 3.46 

AT4G21320 Heat-stress-associated 32-kD protein 4.12 2.18 

AT2G26150 Heat shock transcription factor A2 3.79 3.98 

AT2G25140 HSP100/ClpB4 3.35 3.38 

AT2G19310 HSP20-like chaperones superfamily protein 4.14 2.10 

AT3G12050 Ahal domain-containing protein 2.77 2.25 

AT3G09350 Encodes one of the /Irafa/dops/s Hsp70-binding protein 1 4.12 2.39 

AT2G24860 DnaJ/Hsp40 cysteine-rich domain superfamily protein 6.39 2.12 
Responses to oxidative, salt and other stresses 

AT4G37370 CYP81 D8, cytochrome P450 - like protein 13.24 2.81 

AT2G34500 CYP710A1, putative cytochrome P450 5.73 2.66 

AT1G73480 alpha/beta-Hydrolases superfamily protein 6.47 2.33 

AT3G60140 DIN2, member of glycoside hydrolase family 1 5.49 2.22 

AT5G37710 alpha/beta-Hydrolases superfamily protein 3.33 2.25 

AT3G09640 ATAPX2 5.72 5.17 

AT2G04040 AtDTXI , identified as a detoxifying efflux carrier 4.92 3.00 

AT2G47710 Adenine nucleotide alpha hydrolases-like superfamily 2.35 2.09 
Signal transduction 

AT4G21390 Protein serine/threonine kinase activit/ 17.03 2.00 

AT3G09010 Protein kinase superfamily protein (LRR-RLK VIII) 3.23 2.58 

AT3G22840 Encodes an early light-inducible protein 2.38 2.05 

AT2G21940 Encodes a shikimate kinase 3.62 2.37 

AT2G20900 diacylglycerol kinase 5 (DGK5) 3.57 2.17 
Proteasome related pathway 

AT1G14200 RING/U-box superfamily protein 6.60 2.89 

AT1G26800 RING/U-box superfamily protein 6.22 2.65 

AT1G55530 RING/U-box superfamily protein 5.24 2.45 

AT5G09800 ARM repeat superfamily protein 5.56 2.39 

AT5G48655 RING/U-box superfamily protein 3.95 2.13 
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Table 2 Selection of putative As tolerance-associated genes by comparing the expression ratio between the two 
accessions (Continued) 



AT3G 13430 


RING/U-box superfamily protein 


3.75 


2.28 


AT4G11370 


Encodes a putative RING-H2 finger protein RHAIa 


3.12 


2.15 


AT5G38895 


RING/U-box superfamily protein 


2.83 


2.29 


AT4G23570 


AtSGTIa 


3.55 


2.49 


AT4G 15420 


Ubiquitin fusion degradation UFD1 family protein 


3.43 


2.00 


ATI G32940 


ATSBT3.5, involved in proteolysis 


3.65 


2.05 


AT3G26980 


ML)B4, membrane-anchored ubiquitin-fold protein 4 


4.88 


2.54 


GTP binding 








AT3G21700 


ATSGP2, Monomeric G protein 


4.93 


2.28 


AT1G72660 


GTP1, Small GTP-binding protein 


16.12 


3.01 


Mitochondria electron transport 






AT5G25450 


Cytochrome bd ubiquinol oxidase 


3.91 


2.07 


AT4G28390 


AAC3, ADP.ATP carrier-like protein 


10.16 


2.57 


AT4G27940 


manganese tracking factor for mitochondrial S0D2 (MTMl) 


2.65 


2.70 


Others 








AT2G03430 


Ankyrin repeat family protein 


2.03 


2.37 


AT3G25900 


ATHMT-1, Homocysteine S-methyltransferase 


7.68 


2.04 


AT2G 15490 


ATUGT7, putative glucosyltransferase 


4.86 


2.07 


AT4G00550 


UDP-galactose-dependent digalactosyldiacylglycerol synthase 


6.22 


3.29 


AT4G26270 


phosphofructokinase 3 (ATPFK3) 


5.38 


3.64 


AT3G 16050 


Encodes a protein with pyridoxal phosphate synthase 


3.64 


2.14 


ATI G30070 


SGS domain-containing protein 


3.53 


2.26 


AT5G641 70 


dentin sialophosphoprotein-related protein 


3.29 


2.34 


AT3G57810 


Cysteine proteinases superfamily protein 


2.94 


2.27 


ATI G24090 


RNase H family protein 


2.90 


2.24 


AT2G 16900 


Arabidopsis phospholipase-like protein family 


2.80 


2.08 


AT4G32440 


Plant Tudor-like RNA-binding protein 


2.40 


2.28 


ATI G67360 


Rubber elongation factor protein 


2.24 


2.73 



'Col-0 200/Ws-lOO refers to pair-wise comparison of expression ratio (Col-0 200 \iM As vs Col-0 Control) / (Ws-2 200 As vs Ws-2 Control). 
""Col-O 200/100 refers to pair-wise comparison of expression ratio (Col-O 200 pM As vs Col-0 Control) / (Col-0 100 pM As vs Col-0 Control). 



tolerant and sensitive to As stress. Arabidopsis Col-0 ac- 
cession plants were more tolerant to As stress than Ws- 
2 plants and accumulated significantly less As in root 
tissues. Microarray assay identified the early response of 
Arabidopsis to As stress and possible differences in 
the mechanisms to achieve tolerance in the accessions. 
GO analysis of the transcriptome data suggested that As 
stress significantly affects biological processes related 
to responses to heat, oxidative stress, metal ion and cell 
wall organization, and cytokinin and transport in both 
accessions. The differential expression levels of ABA- 
and ethylene-related genes may contribute to the As 
sensitivity and tolerance of the two accessions. Pathway 
analysis suggested that the As tolerance of the Col-0 
accession was attributed to enhanced expression of regu- 
latory and functional components such as proteasome- 



related mechanisms and responses to heat. The present 
work extends current knowledge of early transcriptional 
regulation by As stress in Arabidopsis roots and provides 
valuable insights into aspects of As toxicity, detoxifica- 
tion and acquired tolerance. 

Ecotypic variation In As tolerance 

Toxicity and tolerance to heavy metals are closely re- 
lated to the accumulation of heavy metals in plant tis- 
sues [9,16,27]. Plant species and even genotypes differ 
greatly in their ability to take up, transport and accumu- 
late heavy metals. In maize, an Al-tolerant genotype 
accumulated significantly less Al in root tips than the 
Al-sensitive genotype [28]. Holcus lanatus has As-tolerant 
and -sensitive accessions, and relatively less As concentra- 
tion accumulated in roots of As-tolerant accessions [23]. 
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Figure 3 Expression analysis of Arabidopsis LRR-RLK VIII gene in 
response to As stress, (a) Semi-quantitative RT-PCR analysis of 
LRR-RLK VIII genes in two Arabidopsis accessions after exposure to 
As stress. The details in annotation of these selected As tolerance- 
associated genes are summarized in Table 2. (b) RT-PCR analysis of 
LRR-RLK VIII genes in response to various stresses from Arabidopsis 
Col-0 accession. The plants were treated with As (100 and 200 ^iM), 
Cu (25 |jM) and H2O2 (100 (iM) for 3, 12 and 24 h. The data are on 
the basis of three biological replicates. 



We selected Arabidopsis with different accessions to study 
As tolerance by root growth inhibition (Figures la and b) 
and As accumulation (Figure Ic). The Arabidopsis acces- 
sion Col-0 was more tolerant to As stress than Ws-2 and 
accumulated less As in root tissues (Figure Ic). An As ex- 
clusion mechanism may operate in roots of As-tolerant 
Col-0. In As-sensitive Ws-2 plants, As toxicity was associ- 
ated with relatively higher accumulation of As in roots 
(Figure 1). GO analysis of microarray data also suggested 
that the increased As accumulation of Ws-2 resulted in 
enhanced expression of toxin and osmotic stress-related 
genes (Table 1). Therefore, the ecotypic variation in As 
tolerance appeared to be associated with As level, as well 
as expression of toxin and osmotic stress-related genes. 

Comparisons of short- and long-term transcrlptome 
responses to As stress 

The transcriptional response of Arabidopsis Col-0 to As (V) 
stress has been reported [10]. As may induce genes involved 
in response to oxidative stress and repress that of genes in- 
duced by phosphate starvation [10]. The cysteine-rich 
metal-binding protein metallothionein (MT) was also in- 
duced. Previous work examined As toxicity in Arabidopsis 



in terms of root growth inhibition and transcriptional re- 
sponses after 3 and 10 d [10]. Forty-six and 113 genes were 
induced and repressed, respectively. The physiological and 
metabolic parameters measured under these long treatment 
periods might be distorted by the severe toxic effects of As. 
Mechanisms of adaptation after long-term exposure are 
relevant. However, we aimed to understand the primary re- 
sponse to metal ion exposure as opposed to responses to 
unspecific cellular damage. Increased transcript abundance 
of heat shock protein and the ubiquitin-proteasome path- 
way (as reported in this study) are general responses to cel- 
lular stress and damage. We examined genomic gene 
expression profiles and biological pathways in Arabidopsis 
in response to short-term As stress (Figures 3 and 4): 1.5 
and 3 h of As exposure resulted in 2451 and 1149 
As-responsive genes in Col-0 and Ws-2, respectively. Only 
16 genes were regulated in common after long-term [10] 
and short-term (this study) As exposure (Additional file 1: 
Table S8). Genes that have previously been reported to be 
As inducible or pathways were not affected by short-term 
As exposure (Table 1). We did not observe changes in ex- 
pression of MT genes, which suggests a unique cellular re- 
sponse to short-term As exposure. Transcrlptome data 
demonstrated that As stress significantly affects biological 
processes related to responses to heat, oxidative stress, 
metal ion, cell wall organization and cytokinin, as well as 
transport (Table 1). Therefore, Arabidopsis plants rapidly 
and simultaneously change the expression of specific sets 
of genes to cope with As stress in root tissues. 

As-responsive genes involved in the detoxification 

Plants possess a range of potential cellular mechanisms 
that may be involved in the detoxification of heavy metals 
such as antioxidant and transporter systems [29,30]. Genes 
coding for proteins involved in oxidative stress such as 
thioredoxins and Class III peroxidase were highly repre- 
sented in our microarray data of Arabidopsis plants ex- 
posed to As stress (Additional file 1: Table S2). Transcripts 
belonging to GST formed the largest group within the 
oxidative stress-related genes. Most of these transcripts 
are in the Tau subfamily (GST-Tau), and the remaining 
sequences are similar to the Phi and lambda subfamily 
(Additional fOe 1: Table S3). GSTs are induced by a num- 
ber of intracellular and environmental factors such as oxi- 
dative stress and heavy metals [31,32] and are involved in 
detoxification of both endogenous and xenobiotic com- 
pounds with electrophilic centers by the nucleophilic 
addition of glutathione [31]. Some isoforms of GST show 
dual activity, additionally functioning as a glutathione per- 
oxidase in the presence of reactive oxygen species [31]. 
This phenomenon provides further evidence for the role 
of GSTs in antioxidant metabolism. As-glutathione conju- 
gates may be produced by GSTs in animal cells such as rat 
liver [15,33,34]. However, GST-mediated conjugation of 
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Figure 4 Functional analysis of a tolerance-associated gene encoding LRR-RLK VIII in response to As stress, (a) Identification of T-DNA 
insertion for LRR-RLK VIII mutations in the Arabldopsis genome by PCR Analysis. The left (LB) border primers for T-DNA insertion, and the left (LP) 
and right (RP) genomic primers for LRR-RLK VIII genes are presented in materials and methods. The sequences of junctions between T-DNA 
borders and the genomic target were detected in homologous LRR-RLK VIII mutants when compared to the wild-type plants. Two T-DNA mutant 
lines for locus AT1G53440 were characterized, and one for AT1G53430 (upper panel). The gene expression of LRR-RLK VIII {AT1G53440) in the 
mutant line (Salk_057812) was analyzed by RT-PCR (bottom panel). There was no characterized T-DNA mutant line for AT3G09010. (b) Effect of As 
(V) on the root elongation of Arabldopsis wild-type and LRR-RLK VIII mutant lines was assessed. Measurement of root elongation was similar to that 
described in Figure 1. Seedlings were grown on quarter-strength MS medium for 4 d and then transferred to medium with 100 and 200 (iM As 
(V) and grown for an additional 4 d. As tolerance was determined by relative root growth after treatment. Root length of plants was measured 
after treatment with As. Root samples were collected from 3 independent experiments (each from a pool of 7 root samples). Data are mean±SD. 
*P < 0.05 compared to As treatment (200 (iM) from wild-type plants. The difference in root elongation is significant according to Student's t test, 
(c) Accumulation of As (V) by Arabldopsis wild-type and LRR-RLK VIII mutant lines was analyzed with the methods similar to Figure 1. Data are 
mean±SD calculated from 3 biological replicates per treatment. *P < 0.05 compared with wild-type plants, (d) Effect of As (200 ^iM), Cu (50 ^iM) 
and H2O2 (100 ^iM) on the root elongation of Arabldopsis wild-type and LRR-RLK VIII mutant lines. Measurement of root elongation was similar to 
that described in Figure 1. At least three biological replicates were performed corresponding to each treatment Data are the mean±SE of three 
independent replicates. 



glutathione with As has not been demonstrated in plants. 
The upregulation of Tau class GSTs has also been noted in 
transcriptomic and proteomic analysis of plant roots under 
As stress [9,32]. Thus, our results suggest that the Tau 
class of GSTs most likely functions at the root level to pro- 
tect cells against heavy-metal — induced oxidative damage. 
In addition, we found differential regulation of genes 
encoding GSTs between As-tolerant and -sensitive Arabl- 
dopsis accessions (Additional file 1: Table S2), with fewer 



As-induced GST genes regulated in As-tolerant Col-0 than 
As-sensitive Ws-2, which indicates the differential GST- 
based detoxification. The GST gene expression profile can 
therefore account for the tolerant Arabldopsis Col-0 acces- 
sion being less affected by As stress. 

As-tolerance associated genes involved in ethylene signalling 

The hormones ethylene participates in signaling cas- 
cades regulating both development and responses to 
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Stress [35]. Recently, cross-talk between ethylene signaling 
and sulfur assimilation in plants has been proposed [36]. 
It was demonstrated that sulfur can induce tolerance to 
Cd stress and alleviate photosynthetic inhibition through 
ethylene by maintaining high GSH levels [37]. However, 
their roles in Arabidopsis exposed to As has not been 
reported. In this study, genes coding for ethylene-signalling 
components were significantly enriched in As-tolerant 
Col-0 with short-term As exposure (Table 1 and Additional 
file 1: Table S2). The transcription factor gene encoding 
ERF (ethylene response factor) was identified as an As- 
tolerance associated gene (AT2G33710) (Table 2). The sul- 
fur assimilation pathway was also affected by As stress 
(Additional file 2: Figure SI). Thus, regulatory interaction 
between sulfur assimilation and ethylene may contribute to 
tolerance to As in Col-0 accession. Furthermore, the pre- 
dominant functional themes of the tolerance-associated 
genes are related to the ubiquitin/proteasome pathway and 
responses to heat (Table 2). These molecular signatures 
identified in the As-responsive transcriptome might work 
in concert for tolerance to As stress. A number of studies 
have suggested a role for ethylene in plants on exposure 
to heavy metals. For example, cadmium or copper in- 
duces the biosynthesis of ethylene in various plant species 
[38,39]. The hormone is transmittable signals, which in- 
voke stress responses and tolerance in the shoot [40,41]. 
In this study, foliar necrosis was less severe in As-tolerant 
Col-0 than As-sensitive Ws-2 after exposure to As for 
1 week (Additional file 4: Figure S3). 

As-tolerance associated genes involved in ubiquitin/ 
proteasome pathways 

The elevated expression of genes related to the proteasome 
pathway indicated that the oxidized or As-bound pro- 
teins in As stressed plant cells might be ubiquitinated for 
proteasome-mediated degradation. Expression levels of the 
components and intermediates of the ubiquitin/proteasome 
pathway were predominately induced in the Col-0 acces- 
sion and identified as tolerant-associated genes under As 
stress (Table 2). A notable number of the proteasome- 
related genes belong to the component of E2 and E3 
ubiquitin ligases (E3s) such as RING/U-box superfarruly 
protein, RING-H2 finger protein and membrane-anchored 
ubiquitin-fold protein (Table 2). In addition, the ubiquitin- 
associated proteins present contained the ubiquitin fusion 
degradation UFDl family protein. These are key intermedi- 
ates in the tagging of proteins for proteasomal degradation. 
RING proteins can act as single components containing 
both the active site and the binding pocket for the E2- 
ubiquitin intermediate [42-44]. Emerging evidence suggests 
the potential roles of the proteasome-mediated pathway 
in plants with As exposure. Expression profiling of Crambe 
abyssinica under As stress identified gene networks in- 
volved in 20S proteasome degradation [27]. In Lemna 



minor, the response to As (III) involves cellular protease- 
chaperone machineries, including heat shock protein 
synthesis and the ubiquitin/proteasome pathway [45]. In 
addition, the As-tolerance associated genes identified in this 
study (Table 2) were compared with Cd-induced genes in 
Arabidopsis [46]. Gene encoding ubiquitin-specific protease 
(AT3G28220) was induced by short-term Cd treatment, in- 
dicating the roles of ubiquitin/proteasome pathways in 
heavy metal stress. Therefore, our results suggest that the 
proteasome-mediated protein degradation machinery is 
vital for As tolerance in plants. 

Potential roles of the As-tolerance associated genes 
encoding LRR-RLKs 

RLKs consist of extracellular repeats that are linked by a 
transmembrane domain to either an intracellular adapter 
domain or a kinase domain [47]. LRR-RLKs represent the 
largest subfamily in the Arabidopsis genome with approxi- 
mately 235 members, divided over 13 subfamilies (LRR I 
to XIII) [48]. The biological functions have been defined 
for only about 30 proteins, which play crucial roles in a 
variety of different physiological processes such as devel- 
opment, pathogen resistance and hormone perception 
[47]. In addition, it has been reported that target genes 
of the Cd-responsive miRNAs encode LRR-RLKs [49]. 
Exploration of Arabidopsis expression data with the Ara- 
bidopsis eFP Browser [50] revealed the LRR-RLK VIII genes 
were upregulated by salt and cold stresses (Additional 
file 5: Figure S4). Functional characterization of the LRR- 
RLK VIII gene (AT1G53430) using a T-DNA insertion 
line showed that the mutant was resistance to salt stress 
[48]. In this study, four genes (AT1G53430, AT1G53440, 
AT3G09010 and AT5G01950) belonging to the LRR-RLK 
VIII family were significantly upregulated by As stress in 
the tolerant Col-0 accession (Table 2 and Additional file 1: 
Table S7). The result was further validated by semi- 
quantitative and real-time RT-PCR analysis (Figure 3). 
Some members of the LRR-RLK VIII family may be func- 
tionally redundant in the regulation of As stress (Figure 4). 
Two of the mutant lines (AT1G53440: SALK_057812 and 
SALK_148231) showed alternation in As sensitivity and 
the metal accumulation (Figure 4). Thus, this study re- 
vealed a novel role for LRR-RLK VIII in sensitivity to As 
stress. Some mammalian toll-like receptors, which are the 
counterparts of plant LRR-RLKs, are involved in xeno- 
biotic response such as drug metabolism [51,52]. It may 
be speculated that specific plant LRR-RLKs could have 
evolved into sensors under heavy metal stress. 

LRR-RLK VIII may involve in perception of As stress 
leading alternation in As sensitivity 

In this study, we report the functional analysis of an 
Arabidopsis LRR-RLK VIII gene involved in the regula- 
tion of early root responses to As stress. The Arabidopsis 
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LRR-RLK VIII mutants alleviated As-induced root growth 
inhibition, indicating that the mutant roots may be less 
sensitive to the heavy metal stress (Figure 4b). Accumula- 
tion of As level was decreased in the LRR-RLK VIII mu- 
tants as compared to wild-type plants (Figure 4c). It is 
suggested that the LRR-RLK VIII mutants are less sensitive 
to As stress due to reduced As accumulation in roots. 
Therefore, the LRR-RLK VIII gene appeared to function as 
a negative regulator of As tolerance. Similarly, knocking 
out expression of a Medicago truncatula LRR-RLK (Srlk) 
by TILLING failed to limit root growth when exposed to 
salt stress [53]. These LRR-RLK mutant plants accumulate 
fewer sodium ions than controls, and several early salt- 
regulated genes are downregulated after exposure to salt 
stress [53]. Taken together, the results may link the LRR- 
RLK receptor with perception of As stress and activa- 
tion of a signaling pathway leading to alternation in As 
sensitivity. 

Conclusions 

We performed a comparative investigation to reveal 
changes in gene expression that take place in roots of an 
As-tolerant and As-sensitive accessions of Arabidopsis 
under short-term As stress. We identified putative candi- 
date genes for As tolerance in Arabidopsis to explore 
molecular mechanisms in response to As. Future studies 
involving Arabidopsis mutants or overexpressors with al- 
tered expression of these genes will be helpful to elucidate 
their biological significance and to clarify new pathways 
involved in As toxicity and tolerance. Analysis of As toler- 
ance by quantitative trait locus (QTL) strategies in succes- 
sive generations obtained from crosses between the two 
accessions wOl facilitate the identification of genes in- 
volved in As tolerance. The knowledge of how plants ac- 
quire tolerance to As is essential for developing novel 
strategies for efficient phytoremediation and As tolerant 
crops. 

Methods 

Plant materials and growth conditions 

Arabidopsis thaliana Col-0 (accession no. CS22625) and 
Ws-2 (no. CS22659) were obtained from the Arabidopsis 
Biological Resource Center (Ohio State University, USA) 
and propagated for 1 generation. The seeds were verna- 
lized for 2 d at 4°C before germination. About 500 seeds 
were sterilised by 10 min of incubation in 1 ml of 1.2% 
(v/v) sodium hypochlorite solution with a few drops of 
Tween 20, then rinsed 4 times with sterilized water. Sur- 
face-sterilized seeds were sown on quarter-strength MS 
(Murashige-Skoog) medium (M0222.0050, Duchefa Bio- 
chemie. The Netherlands) supplemented with 1.0% (w/v) su- 
crose, pH 5.7, with 0.05% (w/v) 2-morpholinoethanesulfonic 
acid monohydrate, and solidified with 1% (w/v) phytagel 
(Sigma, St. Louis, MO, USA). Seedlings were grown under 



16-hr white light/8-hr dark at 25°C with an illumination 
intensity of 3,000 Lux by placing the Petri dish vertically. 

Measurement of root length 

Seven 4-d-old seedlings were transferred to plates (12 x 
12 X 1.7 cm) with medium supplemented with As from 
100 to 300 |iM. Control seedlings were incubated without 
As. Seedlings were allowed to grow for additional 2 d. The 
elongation of primary roots was measured by manually re- 
cording on the plate. Data were obtained from 3 biological 
replicates. 

Determination of As concentration in roots 

After 3 h or 48 h of As treatment, roots of 4-d-old Arabi- 
dopsis seedlings were rinsed thoroughly with distilled 
water and oven dried for 2 d at 60°C, then digested in 
HNO3/H2O2, 3:2 (v/v), for 30 min. As concentration was 
determined by inductively coupled plasma atomic emis- 
sion spectrometry (ICP-AES) (HORIBA Jobin Yvon 2000- 
2, France). Data were mean values from 3 independent 
biological replicates. 

RNA extraction 

The medium was overlaid with a layer of nylon mesh 
and 4-d-old Arabidopsis seedlings grew on the nylon 
mesh to facilitate collection of root samples. Total RNA 
was isolated from roots with use of an RNeasy Plant 
Mini Kit (Qiagen, Hilden, Germany), treated with DNase 
(Roche, Germany) to eliminate genomic DNA contamin- 
ation, then purified and concentrated by use of the 
RNeasy MinElute Cleanup Kit (Qiagen). Three biological 
replicates were grown in the same growth chamber to 
minimize experimental and sample-to-sample variation. 
For each accession, total RNA samples extracted from 
untreated (control) and As-treated (experimental) sam- 
ples were subjected to microarray analysis. 

Microarray analysis 

Microarray analysis involved Affymetrix microarrays 
(GeneChip Arabidopsis ATHl genome array) containing 
22,810 probe sets on a single chip. The data were analyzed 
by use of the Microarray Suite v5.0 with Affjmetrix de- 
fault settings and global scaling as the normalization 
method. The trimmed mean target intensity of each array 
was arbitrarily set to 100. The raw cell intensity data files 
(GeneChip CEL files) were analyzed by use of GeneSpring 
GXIO (AgOent Technologies, CA, USA). Data were nor- 
malized by the GeneChip RMA Robust Multichip Average 
(GC-RMA) algorithm and converted to log2 scale to allow 
for comparing the 3 biological replicates. As treated sample 
was compared to an untreated control sample harvested 
from each accession. Genes with statistically significant 
2-fold change in expression in all experiments were de- 
tected by Student's t test at P < 0.05. The Benjamini and 
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Hochberg algorithm calculates false discovery rates 
(FDRs) that are inherently corrected for multiple testing 
[25]. Genes were considered as significantly up- or 
downregulated if the FDR value for the corresponding 
probe set was < 0.05. We focused our attention on the 
subset of these transcripts that differed by 2- fold be- 
tween As-treated and untreated samples. Three bio- 
logical replicates of each sample were carried out to 
achieve reproducibility of the chip hybridization. The 
microarray data have been deposited in Gene Expres- 
sion Omnibus (no. GEO: GSE31977). 

Gene ontology (GO) analysis 

The differentially expressed genes were classified by bio- 
logical function by analysis with agriGO (http://bioinfo. 
cau.edu.cn/agriGO/) [54] and MapMan [55]. The gene 
lists in each functional categories were obtained from 
The Arabidopsis Information Center (http://www.arabi- 
dopsis.org/) and Arabidopsis thaliana Kinase Database 
(http://bioinformatics.cau.edu.cn/athKD/index.htm). 

Semi-quantitative RT-PCR 

First- strand cDNA was synthesized from 2 |ig of total RNA 
with 1 [il oligo (dT)15 primers by use of the ImProm-II 
Reverse Transcription System (Promega, WI, USA). The 
fragments of cell cycle-related genes were amplified by use 
of gene-specific primers (Additional file 6: Table S9). Actin 
was an internal control. The PGR cycling involved an initial 
denaturation step at 94°C for 2 min, 27-40 cycles of ampli- 
fication and a final elongation step at 72°C for 5 min. PGR 
products were analyzed on a 2% (w/v) agarose gel. The ex- 
periments were repeated at least 2 times for each gene. 

Characterization of Arabidopsis LRR-RLK VIII T-DNA 
insertion mutants 

The Arabidopsis LRR-RLK VIII T-DNA mutant lines 
(AT1G53430, SALK_129312; AT1G53440, SALK_148231 
and SALK_057812; AT3G09010, SALK_019665) were 
obtained from Arabidopsis Biological Resource Genter 
(ABRC; http://abrc.osu.edu/resources). The DNA samples 
were extracted from wild-type and LRR-RLK VIII T-DNA 
mutant lines and subject to PGR analysis. Homozygous 
plants were verified with LRR-RLK V///-specific primers 
spanning the T-DNA insertion site (LP and RP) together 
with a T-DNA left-border primer (LB) (Additional file 6: 
Table S9). The PGR consisted of an initial denaturation 
step at 94°G for 2 min, 30 cycles of amplification with an- 
nealing temperature at 55°G and a final elongation step 
at 72°G for 5 min. The PGR product was separated on a 
2% (w/v) agarose gel. The amplified PGR fragments were 
sequenced to confirm the presence of a LRR-RLK VIII 
sequence and to determine the T-DNA insertion site. 
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